Effects of Bulk Viscosity in Non-linear Bubble Dynamics 
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The non-linear bubble dynamics equations in a compressible liquid have been modified considering 
the effects of compressibility of both the liquid and the gas at the bubble interface. A new bubble 
boundary equation has been derived, which includes a new term resulted from the liquid bulk 
viscosity effects. The influence of this term has been numerically investigated considering the effects 
of water vapor and chemical reactions on the bubble evolution. The results clearly indicate that 
the new term has an important damping role at the collapse, so that its consideration decreases the 
amplitude of the bubble rebounds after the collapse. This damping feature is more remarkable for 
higher deriving pressures. 

PACS numbers: 47.55.Bx, 43.25.Yw, 43.25.+y, 78.60.Mq 



When a small isolated gas bubble, immersed in a liq- 
uid, experiences a high amplitude spherical sound field, 
it grows and contracts non-linearly. Description of the 
dynamics of such non-linear motion is an old challeng- 
ing problem. The radial dynamics of the bubble in an 
incompressible liquid is described by the well-known in- 
compressible Rayleigh-Plesset equation The exten- 
sion of this equation to the bubble motion in a compress- 
ible liquid has been studied by many previous authors 
0, . The most complete existing description was pre- 
sented by Prosperetti and Lezzi [j| . They used a singular- 
perturbation method of the bubble-wall Mach number 
and derived a one-parameter family of equations describ- 
ing the bubble motion in the first order approximation of 
compressibility. This family of equations are written as: 
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where, R, C, Po, P a , and p are bubble radius, liquid 
sound speed, ambient pressure, driving pressure, and 
density of the liquid, respectively. Also, r\ is an arbi- 
trary parameter. Equation Q must be supplemented by 
a boundary condition equation at the bubble interface to 
relate the liquid pressure, Pi, to the gas pressure inside 
the bubble. Like all previous authors, Prosperetti and 
Lezzi used the following incompressible equation for 
this purpose: 
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respectively. Most of the previously obtained equations 
belong to this single parameter family of equations, corre- 
sponding to different values of r\. Moreover, 77 = yields 
results in closest agreement with the numerical simula- 
tion of full partial differential equations 0| . 

In all previous works 0, 0, El 0, , an important ap- 
proximation has been used in the derivation of the bub- 
ble dynamics equations. That is the incompressibility 
assumption of the liquid motion at the bubble interface, 
which has been used in the derivation of Eq'n. (0 ■ Note 
that, all of the effects of the liquid compressibility in all 
previous papers have been resulted from the liquid mo- 
tion around the bubble, but not from the bubble bound- 
ary condition equation. In fact, all previous authors, on 
one hand took into account the compressibility of the 
liquid motion around the bubble, but on the other hand 
neglected its consideration at the bubble interface. 

In this paper, we have modified the bubble dynamics 
equations considering the effects of the liquid compress- 
ibility at the bubble interface. We have derived a new 
bubble boundary equation instead of Eq'n. J2J). The 
new equation has new terms resulted from the effects of 
bulk viscosity of the liquid and the gas. 

To derive the compressible bubble boundary equation, 
the continuity equation and the radial component of the 
stress tensor under the spherical symmetric condition can 
be written as: 
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where, P g , p, and a are gas pressure at the bubble in- 
terface, liquid viscosity coefficient, and surface tension, 



where, p, u, p, and A are density, velocity, pressure, and 
divergence of the velocity, respectively. Also, pb is the 
bulk viscosity coefficient and is defined by Uh = X + 2p/3, 
where A is second coefficient of viscosity 6j. Inserting 
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du/dr from Eq'n. Q, into Eq'n. (@J yields: 

T rr = -p+(/! 6 + ^)A-4^. (5) 
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The velocity divergence, A, can be written as: 

pdt pc 2 dV [ > 

where, the sound speed, c, is defined as c 2 = dp/ dp. The 
boundary continuity requirement at the bubble interface 
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Applying Eq'n. ^ for the gas and the liquid parts of 
Eq'n. JJJ leads to: 
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where, /i g and /Xf, s are the viscosity and the bulk viscos- 
ity coefficients of the gas at the bubble interface, respec- 
tively. Also, A; and A g are the divergence of velocity of 
the liquid and the gas, respectively. Substituting the di- 
vergence of velocity for the liquid and the gas from Eq'n 
© into Eq'n. © yields: 
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where, p g is the gas density at the bubble interface. 
Equation © represents the bubble boundary equation 
containing all effects of the compressibility and the vis- 
cosity of both the liquid and the gas. Comparison of 
Eq'ns. @ and © indicates the existence of three new 
terms in Eq'n. © due to the liquid and the gas com- 
pressibility and viscosity effects. Here, we concentrate 
on the effects of the new term arising from the liquid 
compressibility. Therefore, we neglect the gas viscosity 
because of its smallness relative to the liquid viscosity as 
in previous works Under this circumstance, 

Eq'n. © becomes: 
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It should be mentioned that, although the effects of com- 
pressibility consideration in Eq'n. © are in the first 
order approximation, but these effects have been intro- 
duced completely in Eq'n. Q1U[I. 



To close the mathematical analysis, the gas pressure 
evolution at the bubble interface, P s , must be specified. 
In the most complete approach, it can be determined 
from the simultaneous solution of the conservation equa- 
tions for the bubble interior and the bubble radius equa- 
tions 0>H)E1)E1E1E]- Also, heat conduction and mass 
exchange between the bubble and the surrounding liquid 
affect the bubble evolution. In addition, chemical reac- 
tions occurring in the high temperature conditions at the 
end of the collapse, change the bubble content [l3l fl4|. 
All these complexities have been considered in a complete 
gas dynamics model by Storey and Szeri (lif . 

On the other hand, strong spatial inhomogeneities in- 
side the bubble are not remarkably revealed, unless at 
the end of an intense collapse [ToL HH ]. Therefore, the 
uniformity assumption for the bubble interior seems to 
be useful and provides many features of the bubble mo- 
tion 0, 0| . Using this assumption, recently, Lohse and 
his coworkers presented an ODE model 0, 0, |2(| , in 
which all effects of heat transfer at the bubble interface, 
phase change of water vapor, chemical reactions, and dif- 
fusion of reaction products have been considered. This 
model accurately describes various experimental phase 
diagrams 0] and provides a good agreement with the 
complete direct numerical simulation of Storey and Szeri 

Mm 

Here, for describing the bubble interior evolution, we 
have used the Lohse's group model (the same as what 
has been presented in Ref. |2(j| ) . We do not repeat this 
model here and for more details we refer to Refs. [l9ll20| . 
The calculations were carried out under the framework of 
Eq'n (r) = 0), for both the new compressible (Eq'n. 
I1U[) and the old incompressible (Eq'n. © boundary con- 
ditions. We describe an argon bubble in water at room 
temperature, To = 293.0 K, and atmosphere pressure, 
Pq = 1.0 atm, under the conditions of Single Bubble 
Sonoluminescence |0, The driving pressure was 

P a {t) — P a sin (cot), where u> — 2ir x 26.5 kHz. The con- 
stants and the parameters were set accordingly [ilj ; p = 
998.0 kg/m 3 , C = 1483.0 m/s, fi = 1.01 x 10" 3 kg/ms, 
a = 0.0707 kgs~ 2 . The bulk viscosity of water at room 
temperature was set to be pb = 4.1 x 10 -3 kg/ms [22| . 
The constants and parameters of the gas evolution model 
were set the same as what has been presented in Ref. |20j . 

Figures (1) and (2) illustrate the variations of the 
bubble characteristics (radius, total number of particle 
species, and temperature), for the two boundary condi- 
tion cases. It is observed that the addition of the new vis- 
cous term in Eq'n. I|10|l considerably changes the bubble 
evolution after the collapse. The bubble motion is re- 
markably compressible during the collapse. Therefore, 
the new viscous term, which has been arisen from the 
liquid compressibility, is important in this time interval. 
This term exhibits a damping role and its consideration 
reduces the amplitude of the bubble rebounds. Also, the 
period of the rebounds decreases with the addition of 
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FIG. 1: (a) Time variations of the bubble radius in one period 
according to the new compressible (solid) and the old incom- 
pressible cases, (b) Details of the bubble rebounds after the 
collapse for the two cases, (c) Variation of total number of the 
particle species for the two cases during the bubble rebounds. 
The equilibrium radius is Ro = 3.0 fim and the deriving pres- 
sure is P a — 1.35 atm. 



the new term. Details of our calculations show that the 
minimum radius for the new case is about 10% greater 
than that of the old one. The difference between the two 
cases also appears on the variations of the total number 
of particle species (Ar and H2O plus reactions products) 
after the collapse (Fig. lc). Note that, the difference 
gradually disappears as the bubble rebounds weaken. 

In Fig. (2), details of the gas temperature evolution 
around the minimum radius have been demonstrated. 
Damping feature of the new term is clearly observed by 
a considerable decrease of the peak temperature (about 
50%) and an increase of the temperature pulse width, at 
the collapse time. Also, the time of the peak temperature 
about 5 ns changes with the addition of the new term. 

Figure (3) represents the effects of variations of P a on 
the bubble characteristics at the end of the collapse (peak 
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FIG. 2: Time variations of the gas temperature when the 
bubble reaches its minimum radius according to compress- 
ible (solid) and incompressible (dashed) boundary conditions. 
The parameters and constants are the same as Fig. (1). 



temperature, mole fraction of H2O and reaction prod- 
ucts, and minimum radius), for the two boundary condi- 
tion cases. The ambient radius was fixed (Rq = 5.0 fim). 
Figure 3(a) shows that the peak temperature in the both 
cases grows as the driving pressure is increased. How- 
ever, the rate of increase of the peak temperature for the 
new case is considerably smaller than that of the old one. 
This causes that the difference between the two cases be- 
comes remarkable for the higher driving pressures (about 
50% for P a = 1.5 atm). 

The bubble temperature at the end of the collapse is 
high enough to destroy the chemical bonds of water va- 
por molecules inside the bubble. The products of the dis- 
sociation of water vapor molecules are mainly H2, OH, 
H, O, and O2. The chemical reactions between the par- 
ticle species existing inside the bubble affect the bub- 
ble content at the collapse and its peak temperature 
El El 13 IH El 13 III Here, we have considered 
the effects of the reactions No. l-8ofRefs. The 
dependence of the mole fraction of H2O plus reactions 
products, which is defined as (N tot — Na t ) / Ntot, to the 
driving pressure for the two boundary condition cases has 
been illustrated in Fig. 3(b). It shows that the mole frac- 
tion of H2O plus reactions products is similar for the two 
cases in low amplitudes. The difference appears for the 
higher deriving pressures. It is seen that, for the higher 
deriving pressures, the effect of H2O and reactions prod- 
ucts is more important in the new case relative to the old 
one. 

Figure 3(c) shows the variations of the minimum ra- 
dius as a function of P a , for the two cases. The trend 
of variations is similar for the two cases. But, the min- 
imum radius for the new equation is more than that of 
the old equation because of the reduction of the collapse 
intensity. 

A major deficiency of the old bubble dynamics equa- 
tions is that for strongly driven bubbles, such as sonolu- 
minescence bubbles, large amplitude rebounds are pro- 
duced after the collapse, so that they often last until the 
next acoustic cycle of the periodic driving pressure. This 
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FIG. 3: The bubble characteristics at the time of collapse 
as a function of driving pressure amplitude for the compress- 
ible (solid) and incompressible (dashed) boundary conditions; 
peak temperature (a), mole fraction of H2O and reaction 
products (b), and minimum radius (c). The equilibrium ra- 
dius was fixed (Ro = 5.0 fim). Other constants are the same 
as Figs. (1) and (2). 



is in contrast with the experimental results, which show 
rapidly damped rebounds |2^|. By introducing a damp- 
ing term arisen from the gas compressibility, Moss et. al 
[23j provided a typical solution for this problem. The 
effects of the suggested term by Moss et. al is very simi- 
lar to the damping effects of the new term in this paper, 
(compare Fig. 1(b) with Figs. (3) and (4) of Ref. |23|L 
It seems that the damping feature of the bulk viscosity 
is a better way for solving the mentioned problem. The 
reason is that Eq'n. I|10|) has been derived directly from 
the basic equations of fluid mechanics, on the contrary 
to Eq'n. (3.2) of Ref. |2^|, which was derived by an 
approximate method. 

According to the results of this paper, it is expected 
that the theoretical predictions of the bubble stability 
limits are affected by the addition of the new term to the 



bubble dynamics equations. 
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